In this project we used high‐throughput phenomics to identify synergistic or antagonistic interactions in chemical mixtures in the yeast Saccharomyces cerevisiae and Debaryomyces hansenii and used transcriptomic analysis to understand the mechanistic basis of the interactions.
Transcriptomics analysis to understand the mechanistic basis of chemical interactions was performed for two chemical interactions for growth rate (growth during the exponential growth phase), identified in our phenomics experiment, that were conserved in both species i.e. NaCl and rapamycin synergism, and salt and paraquat antagonism.
The details of the transcriptomic samples are given in COMPILED_DATA/SAMPLE_DETAILS.xlsx
We used the alignment program, HISAT2 (version 2.1.0) Kim, 2015 for mapping next-generation sequencing reads with the reference sequence. Preprocessing of RNA-seq alignments for differential expression calling was performed afterwards using the Python package HTSeq (version 0.11.1) Anders, 2014. The UPSTREAM ANALYSIS RESULTS has been described in details in the Git_folder
Downstream analysis to identify deferentially expressed genes was performed using the R package DESeq2 (version 1.30.1) Love, 2014.
The following packages will be frequently used for the entire data analysis and visualization
#for cerevisiae
Dir_sc <- "RAW_COUNT_DATA/HTSEQ_SC/"
#for hansenii
Dir_dh <- "RAW_COUNT_DATA/HTSEQ_DH/"
#for cerevisiae
sampleinfo_sc <- read.table("COMPILED_DATA/sampleTable_SC.txt", header = TRUE, sep = "\t", as.is = TRUE)
#for hansenii
sampleinfo_dh <- read.table("COMPILED_DATA/sampleTable_DH.txt", header = TRUE, sep = "\t", as.is = TRUE)
#for cerevisiae
sampleTable_sc <- data.frame(sampleName = sampleinfo_sc$sampleName,
fileName = sampleinfo_sc$fileName,
condition = sampleinfo_sc$condition)
#for hansenii
sampleTable_dh <- data.frame(sampleName = sampleinfo_dh$sampleName,
fileName = sampleinfo_dh$fileName,
condition = sampleinfo_dh$condition)
library(DESeq2)
#for cerevisiae
ddsHTSeq_sc <- DESeqDataSetFromHTSeqCount(sampleTable = sampleTable_sc,
directory = Dir_sc,
design= ~ condition)
#for hansenii
ddsHTSeq_dh <- DESeqDataSetFromHTSeqCount(sampleTable = sampleTable_dh,
directory = Dir_dh,
design= ~ condition)
#for cerevisiae
ddsHTSeq_sc
## class: DESeqDataSet
## dim: 7445 18
## metadata(1): version
## assays(1): counts
## rownames(7445): 15S_rRNA 21S_rRNA ... tY(GUA)O tY(GUA)Q
## rowData names(0):
## colnames(18): 1_Sc_C1 2_Sc_C2 ... 17_Sc_SR2 18_Sc_SR3
## colData names(1): condition
#for hansenii
ddsHTSeq_dh
## class: DESeqDataSet
## dim: 6712 21
## metadata(1): version
## assays(1): counts
## rownames(6712): DEHA2A00110g DEHA2A00132g ... ENSRNAG049932294
## ENSRNAG049932295
## rowData names(0):
## colnames(21): 19_Dh_C1 20_Dh_C2 ... 38_Dh_SR2 39_Dh_SR3
## colData names(1): condition
#for cerevisiae
keep<-rowSums(counts(ddsHTSeq_sc)) >= 10
dds_SC<-ddsHTSeq_sc[keep, ]
dds_SC_LOST<-ddsHTSeq_sc[!keep, ]
#for hansenii
keep<-rowSums(counts(ddsHTSeq_dh)) >= 10
dds_DH<-ddsHTSeq_dh[keep, ]
dds_DH_LOST<-ddsHTSeq_dh[!keep, ]
Data-transformation functions within DESeq2 are provided for applications other than differential testing, for example clustering of samples or other machine learning applications. Those are the ntd (shifted logarithm), rlog (regularized logarithm) and vst (variance stabilizing transformation). The shifted logarithm has elevated standard deviation in the lower count range, and the regularized log to a lesser extent, while for the variance stabilized data the standard deviation is roughly constant along the whole dynamic range. Note that the vertical axis in such plots is the square root of the variance over all samples, so including the variance due to the experimental conditions. While a flat curve of the square root of variance over the mean may seem like the goal of such transformations, this may be unreasonable in the case of datasets with many true differences due to the experimental conditions.
#for cerevisiae
rld_sc <- rlogTransformation(dds_SC)
vsd_sc <- varianceStabilizingTransformation(dds_SC)
ntd_sc <- normTransform(dds_SC)
#for hansenii
rld_dh <- rlogTransformation(dds_DH)
vsd_dh <- varianceStabilizingTransformation(dds_DH)
ntd_dh <- normTransform(dds_DH)
Figure 1A: Effects of transformations on the variance, S. cerevisiae and rld transformation
Figure 1B: Effects of transformations on the variance, S. cerevisiae and vsd transformation
Figure 1C: Effects of transformations on the variance, S. cerevisiae and ntd transformation
Figure 2A: Effects of transformations on the variance, D. hansenii and rld transformation
Figure 2B: Effects of transformations on the variance, D. hansenii and vsd transformation
Figure 2C: Effects of transformations on the variance, D. hansenii and ntd transformation
The above analysis shows that the rld data transformation is most suitable for our dataset. Therefore from now on all analysis will be with the rld dataset
Figure 3: PCA plot, rld and S. cerevisiae
Figure 4: PCA plot, rld and D. hansenii
Calculation of the euclidean distance of the transformed data. We need to transpose the data since we want to calculate the distance between samples rather than among genes. So we need the samples as rows rather than columns
library(DESeq2)
library(pheatmap)
library(RColorBrewer)
#for cerevisiae
sampleDists_SC <- dist(t(assay(rld_sc)))
#Creating the corresponding matrix
sampleDistMatrix_SC <- as.matrix(sampleDists_SC)
#creating a distance matrix tagging the groups of the samples
sample_data_sc<-data.frame(Group=sampleinfo_sc$condition)
row.names(sample_data_sc)<-sampleinfo_sc$sampleName
#defining the colors we want to use depending on the condition
annotation_sc <- list(Group = c(Ctrl="green", Para="brown", Rapa="red", Salt="cyan", SaPa="blue", SaRa="magenta"))
#for hansenii
sampleDists_DH <- dist(t(assay(rld_dh)))
#Creating the corresponding matrix
sampleDistMatrix_DH <- as.matrix(sampleDists_DH)
#creating a distance matrix tagging the groups of the samples
sample_data_dh<-data.frame(Group=sampleinfo_dh$condition)
row.names(sample_data_dh)<-sampleinfo_dh$sampleName
#defining the colors we want to use depending on the condition
annotation_dh <- list(Group = c(Ctrl="green", Para="brown", Rapa="red", Salt1="cyan", Salt2="turquoise4", SaPa="blue", SaRa="magenta"))
#Make a color palette for heatmap
colfunc<-colorRampPalette(c("goldenrod4", "goldenrod", "white", "turquoise", "turquoise4"))
Plotting the heatmap for cerevisiae
Figure 5: Similarity matrix, rld and S. cerevisiae
Figure 6: Similarity matrix, rld and D. hansenii
The standard differential expression analysis steps are wrapped into a single function, DESeq.The estimation steps performed by this function are described in the Methods section of the DESeq2 publication Love, 2014.
For differential testing it is recommended to use the DESeq function applied directly to the raw counts.
#for cerevisiae
dds_T_SC <- DESeq(dds_SC)
#for hansenii
dds_T_DH <- DESeq(dds_DH)
#differential expression Salt vs Ctrl
res_SC_Salt <-results(dds_T_SC, contrast=c("condition", "Salt", "Ctrl"))
#differential expression Para vs Ctrl
res_SC_Para <-results(dds_T_SC, contrast=c("condition", "Para", "Ctrl"))
#differential expression SaPa vs Ctrl
res_SC_SaPa <-results(dds_T_SC, contrast=c("condition", "SaPa", "Ctrl"))
#differential expression Rapa vs Ctrl
res_SC_Rapa <-results(dds_T_SC, contrast=c("condition", "Rapa", "Ctrl"))
#differential expression SaRa vs Ctrl
res_SC_SaRa <-results(dds_T_SC, contrast=c("condition", "SaRa", "Ctrl"))
#differential expression SaPa vs Salt
res_SC_MX_SaPavsSalt <-results(dds_T_SC, contrast=c("condition", "SaPa", "Salt"))
#differential expression SaPa vs Para
res_SC_MX_SaPavsPara <-results(dds_T_SC, contrast=c("condition", "SaPa", "Para"))
#differential expression SaRa vs Salt
res_SC_MX_SaRavsSalt <-results(dds_T_SC, contrast=c("condition", "SaRa", "Salt"))
#differential expression SaRa vs Rapa
res_SC_MX_SaRavsRapa <-results(dds_T_SC, contrast=c("condition", "SaRa", "Rapa"))
salt_SC<-data.frame(res_SC_Salt, stringsAsFactors = FALSE)
para_SC<-data.frame(res_SC_Para, stringsAsFactors = FALSE)
rapa_SC<-data.frame(res_SC_Rapa, stringsAsFactors = FALSE)
SaPa_SC<-data.frame(res_SC_SaPa, stringsAsFactors = FALSE)
SaRa_SC<-data.frame(res_SC_SaRa, stringsAsFactors = FALSE)
SaPavsSalt_SC <- data.frame(res_SC_MX_SaPavsSalt, stringsAsFactors = FALSE)
SaPavsPara_SC <- data.frame(res_SC_MX_SaPavsPara, stringsAsFactors = FALSE)
SaRavsSalt_SC <- data.frame(res_SC_MX_SaRavsSalt, stringsAsFactors = FALSE)
SaRavsRapa_SC <- data.frame(res_SC_MX_SaRavsRapa, stringsAsFactors = FALSE)
#Import the Gene description file named **Gene_List_SC.tsv** in the **COMPILED_DATA** folder
Gene_Description_SC <- read.csv(file = "COMPILED_DATA/Gene_List_SC.tsv",
header = FALSE,
stringsAsFactors = FALSE,
sep = "\t", quote = "", fill = FALSE)
colnames(Gene_Description_SC) <- c("SGD_DB_ID",
"SYS_ID",
"ORG",
"GENE_SYM",
"NAME",
"DESCRIPTION")
row.names(Gene_Description_SC) <- Gene_Description_SC$SYS_ID
#Connect the result data.frame with the gene description
y <- as.character(rownames(salt_SC))
x <- Gene_Description_SC[(Gene_Description_SC$SYS_ID %in% y), ]
#column binding with the dataset
salt_SC_des <- cbind(salt_SC[x$SYS_ID, ], x)
para_SC_des <- cbind(para_SC[x$SYS_ID, ], x)
rapa_SC_des <- cbind(rapa_SC[x$SYS_ID, ], x)
SaPa_SC_des <- cbind(SaPa_SC[x$SYS_ID, ], x)
SaRa_SC_des <- cbind(SaRa_SC[x$SYS_ID, ], x)
SaPavsSalt_SC_des <- cbind(SaPavsSalt_SC[x$SYS_ID, ], x)
SaPavsPara_SC_des <- cbind(SaPavsPara_SC[x$SYS_ID, ], x)
SaRavsSalt_SC_des <- cbind(SaRavsSalt_SC[x$SYS_ID, ], x)
SaRavsRapa_SC_des <- cbind(SaRavsRapa_SC[x$SYS_ID, ], x)
#Write table E.g.
#write.table(SaRa_SC_des, file = "COMPILED_DATA/SARAvsCTRL_SC.tsv", sep = "\t", na = "NA", row.names = FALSE, col.names = TRUE)
#differential expression Salt1 vs Ctrl
res_DH_Salt1 <-results(dds_T_DH, contrast=c("condition", "Salt1", "Ctrl"))
#differential expression Salt2 vs Ctrl
res_DH_Salt2 <-results(dds_T_DH, contrast=c("condition", "Salt2", "Ctrl"))
#differential expression Para vs Ctrl
res_DH_Para <-results(dds_T_DH, contrast=c("condition", "Para", "Ctrl"))
#differential expression SaPa vs Ctrl
res_DH_SaPa <-results(dds_T_DH, contrast=c("condition", "SaPa", "Ctrl"))
#differential expression Rapa vs Ctrl
res_DH_Rapa <-results(dds_T_DH, contrast=c("condition", "Rapa", "Ctrl"))
#differential expression SaRa vs Ctrl
res_DH_SaRa <-results(dds_T_DH, contrast=c("condition", "SaRa", "Ctrl"))
#differential expression SaPa vs Salt1
res_DH_MX_SaPavsSalt1 <-results(dds_T_DH, contrast=c("condition", "SaPa", "Salt1"))
#differential expression SaPa vs Para
res_DH_MX_SaPavsPara <-results(dds_T_DH, contrast=c("condition", "SaPa", "Para"))
#differential expression SaRa vs Salt
res_DH_MX_SaRavsSalt2 <-results(dds_T_DH, contrast=c("condition", "SaRa", "Salt2"))
#differential expression SaRa vs Rapa
res_DH_MX_SaRavsRapa <-results(dds_T_DH, contrast=c("condition", "SaRa", "Rapa"))
salt1_DH<-data.frame(res_DH_Salt1, stringsAsFactors = FALSE)
salt2_DH<-data.frame(res_DH_Salt2, stringsAsFactors = FALSE)
para_DH<-data.frame(res_DH_Para, stringsAsFactors = FALSE)
rapa_DH<-data.frame(res_DH_Rapa, stringsAsFactors = FALSE)
SaPa_DH<-data.frame(res_DH_SaPa, stringsAsFactors = FALSE)
SaRa_DH<-data.frame(res_DH_SaRa, stringsAsFactors = FALSE)
SaPavsSalt1_DH <- data.frame(res_DH_MX_SaPavsSalt1, stringsAsFactors = FALSE)
SaPavsPara_DH <- data.frame(res_DH_MX_SaPavsPara, stringsAsFactors = FALSE)
SaRavsSalt2_DH <- data.frame(res_DH_MX_SaRavsSalt2, stringsAsFactors = FALSE)
SaRavsRapa_DH <- data.frame(res_DH_MX_SaRavsRapa, stringsAsFactors = FALSE)
Gene decsription and GO annotations for D. hansenii was downloaded from uniprot
Orthologous gene analysis between S. cerevisiae and D. hansenii was performed using Orthovenn2, Ling Xu et al. 2019
The orthologs output from the Orthovenn2 was manually curated in the following way to connect it with the Gene description list obtained from uniprot.
#Import the Gene description file named **uniprot_DH_CBS767_Gene_List.tab** in the **COMPILED_DATA** folder
Gene_Description_DH <- read.csv(file = "COMPILED_DATA/uniprot_DH_CBS767_Gene_List.tab",
header = TRUE,
stringsAsFactors = FALSE,
sep = "\t", quote = "", fill = FALSE)
colnames(Gene_Description_DH) <- c("UNIPROT_ENTRY",
"STATUS",
"PROTEIN_NAME",
"ORG",
"PROTEIN_LENGTH",
"SYS_ID",
"GENE_SYM",
"KEGG_CR",
"GO_ALL",
"GO_PROCESS",
"GO_COMPONENT",
"GO_FUNCTION")
Gene_Description_DH <- Gene_Description_DH[-which(duplicated(Gene_Description_DH$SYS_ID)), ]
row.names(Gene_Description_DH) <- Gene_Description_DH$SYS_ID
# Adding the orthologs of S. cerevisiae to this data.frame
# We Extracted the Orthologous gene cluster analysis between S. cerevisiae and D. hansenii using Orthovenn2
Ortho <- read.table("COMPILED_DATA/Ortholog_DHvsSC.txt", header = FALSE, sep = "\t", as.is = TRUE)
colnames(Ortho) <- c("DH_SYS_ID", "SC_SYS_ID", "ORTHO_SCORE")
x <- data.frame()
for(i in 1:nrow(Gene_Description_DH)){
x <- Ortho[which(Ortho$DH_SYS_ID==Gene_Description_DH$SYS_ID[i]), ]
if(nrow(x)==0){
Gene_Description_DH[i, 13:14] <- NA
} else{
Gene_Description_DH[i, 13:14] <- x[1, 2:3]
}
}
colnames(Gene_Description_DH)[13:14] <- c("SC_SYS_ID", "ORTHO_SCORE")
# Next, we connect some information of the SC ortholog of DH gene
x <- data.frame()
for(i in 1:nrow(Gene_Description_DH)){
x <- Gene_Description_SC[which(Gene_Description_SC$SYS_ID==Gene_Description_DH$SC_SYS_ID[i]), ]
if(nrow(x)==0){
Gene_Description_DH[i, 15:18] <- NA
} else{
Gene_Description_DH[i, 15:18] <- x[1, c(1, 4:6)]
}
}
colnames(Gene_Description_DH)[15:18] <- c("SGD_DB_ID", "SC_GENE_SYM", "NAME", "DESCRIPTION")
#Connecting the dataset
y <- as.character(rownames(salt1_DH))
x <- Gene_Description_DH[(Gene_Description_DH$SYS_ID %in% y), ]
#column binding with the dataset
salt1_DH_des <- cbind(salt1_DH[x$SYS_ID, ], x)
salt2_DH_des <- cbind(salt2_DH[x$SYS_ID, ], x)
para_DH_des <- cbind(para_DH[x$SYS_ID, ], x)
rapa_DH_des <- cbind(rapa_DH[x$SYS_ID, ], x)
SaPa_DH_des <- cbind(SaPa_DH[x$SYS_ID, ], x)
SaRa_DH_des <- cbind(SaRa_DH[x$SYS_ID, ], x)
SaPavsSalt1_DH_des <- cbind(SaPavsSalt1_DH[x$SYS_ID, ], x)
SaPavsPara_DH_des <- cbind(SaPavsPara_DH[x$SYS_ID, ], x)
SaRavsSalt2_DH_des <- cbind(SaRavsSalt2_DH[x$SYS_ID, ], x)
SaRavsRapa_DH_des <- cbind(SaRavsRapa_DH[x$SYS_ID, ], x)
#Write table E.g.
#write.table(salt1_DH_des, file = "COMPILED_DATA/DEG_DATA_DH/SALT1vsCTRL_DH.tsv", sep = "\t", na = "NA", row.names = FALSE, col.names = TRUE)
Figure 7: Histogram, Single compound effect size
A scatter plot of log2 fold changes (on the y-axis) versus the mean of normalized counts (on the x-axis)
Shrinkage of effect size (LFC estimates) is useful for visualization and ranking of genes. To shrink the LFC, we pass the dds object to the function lfcShrink. Below we specify to use the apeglm method for effect size shrinkage Zhu, Ibrahim, and Love 2018, which improves on the previous estimator.
We provide the dds object and the name or number of the coefficient we want to shrink, where the number refers to the order of the coefficient as it appears after using the resultsNames() function
#
resultsNames(dds_T_SC)
## [1] "Intercept" "condition_Para_vs_Ctrl" "condition_Rapa_vs_Ctrl"
## [4] "condition_Salt_vs_Ctrl" "condition_SaPa_vs_Ctrl" "condition_SaRa_vs_Ctrl"
resultsNames(dds_T_DH)
## [1] "Intercept" "condition_Para_vs_Ctrl"
## [3] "condition_Rapa_vs_Ctrl" "condition_Salt1_vs_Ctrl"
## [5] "condition_Salt2_vs_Ctrl" "condition_SaPa_vs_Ctrl"
## [7] "condition_SaRa_vs_Ctrl"
#Data_preparation_SC
resLFC_SC_Salt <- lfcShrink(dds_T_SC, coef = "condition_Salt_vs_Ctrl")
resLFC_SC_Para <- lfcShrink(dds_T_SC, coef = "condition_Para_vs_Ctrl")
resLFC_SC_SaPa <- lfcShrink(dds_T_SC, coef = "condition_SaPa_vs_Ctrl")
resLFC_SC_Rapa <- lfcShrink(dds_T_SC, coef = "condition_Rapa_vs_Ctrl")
resLFC_SC_SaRa <- lfcShrink(dds_T_SC, coef = "condition_SaRa_vs_Ctrl")
#Data_preparation_DH
resLFC_DH_Salt1 <- lfcShrink(dds_T_DH, coef = "condition_Salt1_vs_Ctrl")
resLFC_DH_Salt2 <- lfcShrink(dds_T_DH, coef = "condition_Salt2_vs_Ctrl")
resLFC_DH_Para <- lfcShrink(dds_T_DH, coef = "condition_Para_vs_Ctrl")
resLFC_DH_SaPa <- lfcShrink(dds_T_DH, coef = "condition_SaPa_vs_Ctrl")
resLFC_DH_Rapa <- lfcShrink(dds_T_DH, coef = "condition_Rapa_vs_Ctrl")
resLFC_DH_SaRa <- lfcShrink(dds_T_DH, coef = "condition_SaRa_vs_Ctrl")
A. UNSHRUNKEN
Figure 8: plotMA unshrunken, Single compound effect size
B. SHRUNKEN USING APEGLM
Figure 9: plotMA shrunken, Single compound effect size
A. S. CEREVISIAE SINGLE COMPOUND
## Loading required package: ggplot2
## Loading required package: ggrepel
## Registered S3 methods overwritten by 'ggalt':
## method from
## grid.draw.absoluteGrob ggplot2
## grobHeight.absoluteGrob ggplot2
## grobWidth.absoluteGrob ggplot2
## grobX.absoluteGrob ggplot2
## grobY.absoluteGrob ggplot2
## Warning: One or more p-values is 0. Converting to 10^-1 * current lowest non-
## zero p-value...
Figure 10: Volcano plot- SALT vs CTRL, S. cerevisiae, Single compound effect size
## Warning: One or more p-values is 0. Converting to 10^-1 * current lowest non-
## zero p-value...
Figure 11: Volcano plot- PARA vs CTRL, S. cerevisiae, Single compound effect size
## Warning: One or more p-values is 0. Converting to 10^-1 * current lowest non-
## zero p-value...
Figure 12: Volcano plot- RAPA vs CTRL, S. cerevisiae, Single compound effect size
We performed GO enrichment analysis for both species for
# WITH GROWTH RELATED GENES
GO_salt_SC_all <- salt_SC_des[which(abs(salt_SC_des$log2FoldChange) >= 1
& salt_SC_des$padj <= 0.05), 7]
GO_salt_SC_up <- salt_SC_des[which(salt_SC_des$log2FoldChange>= 1
& salt_SC_des$padj <= 0.05), 7]
GO_salt_SC_down <- salt_SC_des[which(salt_SC_des$log2FoldChange <= -1
& salt_SC_des$padj <= 0.05), 7]
# WITHOUT GROWTH RELATED GENES
#Extracting GROWTH RELATED GENES SGD ID
GRG_ALL_SGD_DB_ID <- salt_SC_des$SGD_DB_ID[which(salt_SC_des$SYS_ID %in% c(Gene_names_ex_up, Gene_names_ex_down))]
# Writing the Background set of SGD_ID i.e. all SGD_ID without Growth related gene IDs
BACKGROUND_GENESET_SGD_DB_ID <- salt_SC_des[which(!(salt_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
# Extracting list for GO analysis
GO_salt_SC_all_wo_GRG <- salt_SC_des[which((abs(salt_SC_des$log2FoldChange) >= 1 & salt_SC_des$padj <= 0.05)
& !(salt_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_salt_SC_up_wo_GRG <- salt_SC_des[which((salt_SC_des$log2FoldChange >= 1 & salt_SC_des$padj <= 0.05)
& !(salt_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_salt_SC_down_wo_GRG <- salt_SC_des[which((salt_SC_des$log2FoldChange <= -1 & salt_SC_des$padj <= 0.05)
& !(salt_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
#WRITE TABLE
#Eg. write.table(GO_salt_SC_all, file = "COMPILED_DATA/GO_GENE_LISTS/GO_salt_SC_all.txt", sep = "\t", na = "NA", row.names = FALSE, col.names = FALSE, quote = FALSE)
#write.table(BACKGROUND_GENESET_SGD_DB_ID, file = "COMPILED_DATA/GO_GENE_LISTS/BACKGROUND_GENESET_SGD_DB_ID.txt", sep = "\t", na = "NA", row.names = FALSE, col.names = FALSE, quote = FALSE)
# WITH GROWTH RELATED GENES
GO_para_SC_all <- para_SC_des[which(abs(para_SC_des$log2FoldChange) >= 1
& para_SC_des$padj <= 0.05), 7]
GO_para_SC_up <- para_SC_des[which(para_SC_des$log2FoldChange>= 1
& para_SC_des$padj <= 0.05), 7]
GO_para_SC_down <- para_SC_des[which(para_SC_des$log2FoldChange <= -1
& para_SC_des$padj <= 0.05), 7]
# WITHOUT GROWTH RELATED GENES
GO_para_SC_all_wo_GRG <- para_SC_des[which((abs(para_SC_des$log2FoldChange) >= 1 & para_SC_des$padj <= 0.05)
& !(para_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_para_SC_up_wo_GRG <- para_SC_des[which((para_SC_des$log2FoldChange >= 1 & para_SC_des$padj <= 0.05)
& !(para_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_para_SC_down_wo_GRG <- para_SC_des[which((para_SC_des$log2FoldChange <= -1 & para_SC_des$padj <= 0.05)
& !(para_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
#WRITE TABLE
#Eg. write.table(GO_para_SC_all, file = "COMPILED_DATA/GO_GENE_LISTS/GO_para_SC_all.txt", sep = "\t", na = "NA", row.names = FALSE, col.names = FALSE, quote = FALSE)
# WITH GROWTH RELATED GENES
GO_rapa_SC_all <- rapa_SC_des[which(abs(rapa_SC_des$log2FoldChange) >= 1
& rapa_SC_des$padj <= 0.05), 7]
GO_rapa_SC_up <- rapa_SC_des[which(rapa_SC_des$log2FoldChange>= 1
& rapa_SC_des$padj <= 0.05), 7]
GO_rapa_SC_down <- rapa_SC_des[which(rapa_SC_des$log2FoldChange <= -1
& rapa_SC_des$padj <= 0.05), 7]
# WITHOUT GROWTH RELATED GENES
GO_rapa_SC_all_wo_GRG <- rapa_SC_des[which((abs(rapa_SC_des$log2FoldChange) >= 1 & rapa_SC_des$padj <= 0.05)
& !(rapa_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_rapa_SC_up_wo_GRG <- rapa_SC_des[which((rapa_SC_des$log2FoldChange >= 1 & rapa_SC_des$padj <= 0.05)
& !(rapa_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_rapa_SC_down_wo_GRG <- rapa_SC_des[which((rapa_SC_des$log2FoldChange <= -1 & rapa_SC_des$padj <= 0.05)
& !(rapa_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
#WRITE TABLE
#Eg. write.table(GO_rapa_SC_all, file = "COMPILED_DATA/GO_GENE_LISTS/GO_rapa_SC_all.txt", sep = "\t", na = "NA", row.names = FALSE, col.names = FALSE, quote = FALSE)
# WITH GROWTH RELATED GENES
GO_SaPa_SC_all <- SaPa_SC_des[which(abs(SaPa_SC_des$log2FoldChange) >= 1
& SaPa_SC_des$padj <= 0.05), 7]
GO_SaPa_SC_up <- SaPa_SC_des[which(SaPa_SC_des$log2FoldChange>= 1
& SaPa_SC_des$padj <= 0.05), 7]
GO_SaPa_SC_down <- SaPa_SC_des[which(SaPa_SC_des$log2FoldChange <= -1
& SaPa_SC_des$padj <= 0.05), 7]
# WITHOUT GROWTH RELATED GENES
GO_SaPa_SC_all_wo_GRG <- SaPa_SC_des[which((abs(SaPa_SC_des$log2FoldChange) >= 1 & SaPa_SC_des$padj <= 0.05)
& !(SaPa_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_SaPa_SC_up_wo_GRG <- SaPa_SC_des[which((SaPa_SC_des$log2FoldChange >= 1 & SaPa_SC_des$padj <= 0.05)
& !(SaPa_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_SaPa_SC_down_wo_GRG <- SaPa_SC_des[which((SaPa_SC_des$log2FoldChange <= -1 & SaPa_SC_des$padj <= 0.05)
& !(SaPa_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
#WRITE TABLE
#Eg. write.table(GO_SaPa_SC_all, file = "COMPILED_DATA/GO_GENE_LISTS/GO_SaPa_SC_all.txt", sep = "\t", na = "NA", row.names = FALSE, col.names = FALSE, quote = FALSE)
# WITH GROWTH RELATED GENES
GO_SaRa_SC_all <- SaRa_SC_des[which(abs(SaRa_SC_des$log2FoldChange) >= 1
& SaRa_SC_des$padj <= 0.05), 7]
GO_SaRa_SC_up <- SaRa_SC_des[which(SaRa_SC_des$log2FoldChange>= 1
& SaRa_SC_des$padj <= 0.05), 7]
GO_SaRa_SC_down <- SaRa_SC_des[which(SaRa_SC_des$log2FoldChange <= -1
& SaRa_SC_des$padj <= 0.05), 7]
# WITHOUT GROWTH RELATED GENES
GO_SaRa_SC_all_wo_GRG <- SaRa_SC_des[which((abs(SaRa_SC_des$log2FoldChange) >= 1 & SaRa_SC_des$padj <= 0.05)
& !(SaRa_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_SaRa_SC_up_wo_GRG <- SaRa_SC_des[which((SaRa_SC_des$log2FoldChange >= 1 & SaRa_SC_des$padj <= 0.05)
& !(SaRa_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_SaRa_SC_down_wo_GRG <- SaRa_SC_des[which((SaRa_SC_des$log2FoldChange <= -1 & SaRa_SC_des$padj <= 0.05)
& !(SaRa_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
#WRITE TABLE
#Eg. write.table(GO_SaRa_SC_all, file = "COMPILED_DATA/GO_GENE_LISTS/GO_SaRa_SC_all.txt", sep = "\t", na = "NA", row.names = FALSE, col.names = FALSE, quote = FALSE)
# WITH GROWTH RELATED GENES
GO_SaPavsSalt_SC_all <- SaPavsSalt_SC_des[which(abs(SaPavsSalt_SC_des$log2FoldChange) >= 1
& SaPavsSalt_SC_des$padj <= 0.05), 7]
GO_SaPavsSalt_SC_up <- SaPavsSalt_SC_des[which(SaPavsSalt_SC_des$log2FoldChange>= 1
& SaPavsSalt_SC_des$padj <= 0.05), 7]
GO_SaPavsSalt_SC_down <- SaPavsSalt_SC_des[which(SaPavsSalt_SC_des$log2FoldChange <= -1
& SaPavsSalt_SC_des$padj <= 0.05), 7]
# WITHOUT GROWTH RELATED GENES
GO_SaPavsSalt_SC_all_wo_GRG <- SaPavsSalt_SC_des[which((abs(SaPavsSalt_SC_des$log2FoldChange) >= 1 & SaPavsSalt_SC_des$padj <= 0.05)
& !(SaPavsSalt_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_SaPavsSalt_SC_up_wo_GRG <- SaPavsSalt_SC_des[which((SaPavsSalt_SC_des$log2FoldChange >= 1 & SaPavsSalt_SC_des$padj <= 0.05)
& !(SaPavsSalt_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_SaPavsSalt_SC_down_wo_GRG <- SaPavsSalt_SC_des[which((SaPavsSalt_SC_des$log2FoldChange <= -1 & SaPavsSalt_SC_des$padj <= 0.05)
& !(SaPavsSalt_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
#WRITE TABLE
#Eg. write.table(GO_SaPavsSalt_SC_all, file = "COMPILED_DATA/GO_GENE_LISTS/GO_SaPavsSalt_SC_all.txt", sep = "\t", na = "NA", row.names = FALSE, col.names = FALSE, quote = FALSE)
# WITH GROWTH RELATED GENES
GO_SaPavsPara_SC_all <- SaPavsPara_SC_des[which(abs(SaPavsPara_SC_des$log2FoldChange) >= 1
& SaPavsPara_SC_des$padj <= 0.05), 7]
GO_SaPavsPara_SC_up <- SaPavsPara_SC_des[which(SaPavsPara_SC_des$log2FoldChange>= 1
& SaPavsPara_SC_des$padj <= 0.05), 7]
GO_SaPavsPara_SC_down <- SaPavsPara_SC_des[which(SaPavsPara_SC_des$log2FoldChange <= -1
& SaPavsPara_SC_des$padj <= 0.05), 7]
# WITHOUT GROWTH RELATED GENES
GO_SaPavsPara_SC_all_wo_GRG <- SaPavsPara_SC_des[which((abs(SaPavsPara_SC_des$log2FoldChange) >= 1 & SaPavsPara_SC_des$padj <= 0.05)
& !(SaPavsPara_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_SaPavsPara_SC_up_wo_GRG <- SaPavsPara_SC_des[which((SaPavsPara_SC_des$log2FoldChange >= 1 & SaPavsPara_SC_des$padj <= 0.05)
& !(SaPavsPara_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_SaPavsPara_SC_down_wo_GRG <- SaPavsPara_SC_des[which((SaPavsPara_SC_des$log2FoldChange <= -1 & SaPavsPara_SC_des$padj <= 0.05)
& !(SaPavsPara_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
#WRITE TABLE
#Eg. write.table(GO_SaPavsPara_SC_all, file = "COMPILED_DATA/GO_GENE_LISTS/GO_SaPavsPara_SC_all.txt", sep = "\t", na = "NA", row.names = FALSE, col.names = FALSE, quote = FALSE)
# WITH GROWTH RELATED GENES
GO_SaRavsSalt_SC_all <- SaRavsSalt_SC_des[which(abs(SaRavsSalt_SC_des$log2FoldChange) >= 1
& SaRavsSalt_SC_des$padj <= 0.05), 7]
GO_SaRavsSalt_SC_up <- SaRavsSalt_SC_des[which(SaRavsSalt_SC_des$log2FoldChange>= 1
& SaRavsSalt_SC_des$padj <= 0.05), 7]
GO_SaRavsSalt_SC_down <- SaRavsSalt_SC_des[which(SaRavsSalt_SC_des$log2FoldChange <= -1
& SaRavsSalt_SC_des$padj <= 0.05), 7]
# WITHOUT GROWTH RELATED GENES
GO_SaRavsSalt_SC_all_wo_GRG <- SaRavsSalt_SC_des[which((abs(SaRavsSalt_SC_des$log2FoldChange) >= 1 & SaRavsSalt_SC_des$padj <= 0.05)
& !(SaRavsSalt_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_SaRavsSalt_SC_up_wo_GRG <- SaRavsSalt_SC_des[which((SaRavsSalt_SC_des$log2FoldChange >= 1 & SaRavsSalt_SC_des$padj <= 0.05)
& !(SaRavsSalt_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_SaRavsSalt_SC_down_wo_GRG <- SaRavsSalt_SC_des[which((SaRavsSalt_SC_des$log2FoldChange <= -1 & SaRavsSalt_SC_des$padj <= 0.05)
& !(SaRavsSalt_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
#WRITE TABLE
#Eg. write.table(GO_SaRavsSalt_SC_all, file = "COMPILED_DATA/GO_GENE_LISTS/GO_SaRavsSalt_SC_all.txt", sep = "\t", na = "NA", row.names = FALSE, col.names = FALSE, quote = FALSE)
# WITH GROWTH RELATED GENES
GO_SaRavsRapa_SC_all <- SaRavsRapa_SC_des[which(abs(SaRavsRapa_SC_des$log2FoldChange) >= 1
& SaRavsRapa_SC_des$padj <= 0.05), 7]
GO_SaRavsRapa_SC_up <- SaRavsRapa_SC_des[which(SaRavsRapa_SC_des$log2FoldChange>= 1
& SaRavsRapa_SC_des$padj <= 0.05), 7]
GO_SaRavsRapa_SC_down <- SaRavsRapa_SC_des[which(SaRavsRapa_SC_des$log2FoldChange <= -1
& SaRavsRapa_SC_des$padj <= 0.05), 7]
# WITHOUT GROWTH RELATED GENES
GO_SaRavsRapa_SC_all_wo_GRG <- SaRavsRapa_SC_des[which((abs(SaRavsRapa_SC_des$log2FoldChange) >= 1 & SaRavsRapa_SC_des$padj <= 0.05)
& !(SaRavsRapa_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_SaRavsRapa_SC_up_wo_GRG <- SaRavsRapa_SC_des[which((SaRavsRapa_SC_des$log2FoldChange >= 1 & SaRavsRapa_SC_des$padj <= 0.05)
& !(SaRavsRapa_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
GO_SaRavsRapa_SC_down_wo_GRG <- SaRavsRapa_SC_des[which((SaRavsRapa_SC_des$log2FoldChange <= -1 & SaRavsRapa_SC_des$padj <= 0.05)
& !(SaRavsRapa_SC_des$SGD_DB_ID %in% GRG_ALL_SGD_DB_ID)), 7]
#WRITE TABLE
#Eg. write.table(GO_SaRavsRapa_SC_all, file = "COMPILED_DATA/GO_GENE_LISTS/GO_SaRavsRapa_SC_all.txt", sep = "\t", na = "NA", row.names = FALSE, col.names = FALSE, quote = FALSE)